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Abstract 

From the 16th to the 18th centuries in Japan, saltpeter was produced using a biological niter-bed process and was formed 
under the floor of gassho-style houses in the historic villages of Shirakawa-go and Gokayama, which are classified as United 
Nations Educational, Scientific and Cultural Organization (UNESCO) World Heritage Sites. The relict niter-beds are now 
conserved in the underfloor space of gassho-style houses, where they are isolated from destabilizing environmental factors 
and retain the ability to produce nitrate. However, little is known about the nitrifying microbes in such relict niter-bed 
ecosystems. In this study, the microbial community structures within nine relict niter-bed soils were investigated using 454 
pyrotag analysis targeting the 16S rRNA gene and the bacterial and archaeal ammonia monooxygenase gene (amoA). The 
16S rRNA gene pyrotag analysis showed that members of the phyla Proteobacteria, Actinobacteria, Bacteroidetes, Chloroflexi, 
Firmicutes, Gemmatimonadetes, and Planctomycetes were major microbial constituents, and principal coordinate analysis 
showed that the N0 3 ~, Cl~, K + , and Na + contents were potential determinants of the structures of entire microbial 
communities in relict niter-bed soils. The bacterial and archaeal amoA libraries indicated that members of the Nitrosospira- 
type ammonia-oxidizing bacteria (AOB) and "Ca. Nitrososphaera"-type ammonia-oxidizing archaea (AOA), respectively, 
predominated in relict niter-bed soils. In addition, soil pH and organic carbon content were important factors for the 
ecological niche of AOB and AOA in relict niter-bed soil ecosystems. 
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Introduction 

Until the early 20th century, saltpeter (i.e., niter or potassium 
nitrate [KNO :i ]), which is known as the main source of 
gunpowder, was produced via an artificial bioprocess, the so- 
called "niter-bed" [1,2]. The basic concept of saltpeter production 
via niter-bed involves the nitrification of ammonia-containing 
wastes (e.g., plant residue, horse manure, and human urine) mixed 
with a certain amount of soil [1,2]. This nitrification reaction was 
thought to be performed by ammonia-oxidizing and nitrite- 
oxidizing microbes [1,2]. The first description of a niter-bed was 
found in the military technological manual "Bellifortis," which was 
written by Konrad Kyeser and published in 1405. This manual 
contains brief instructions describing how to construct a niter-bed 
[3,4], and this process was employed for saltpeter production in 
European countries [1,5,6], India [7,8], and Japan [9,10] until the 
discovery of a large natural deposit in Chile [11]. In Finland, the 
export of saltpeter produced by niter-bed reached a maximum 
level (433 tons) in 1816 during the last stages of the Napoleonic 
wars [6]. In India, a large amount of saltpeter was produced 



annually (20,000 tons/year) using the niter-bed methodology in 
the early 20th century [7]. Thus, the niter-bed process is a 
historically important industry in several nations. 

From the 16th to the 18th centuries in Japan, saltpeter was 
produced for military purposes using the niter-bed formed under 
the floor of gassho-style houses [10], which is a specific 
architectural style used in the historic villages of Shirakawa-go 
and Gokayama, which are listed as UNESCO World Heritage 
Sites [12]. Because these villages were located in the mountainous 
area of the former Kaga domain in Japan (currently the Toyama 
and Ishikawa prefectures), there was limited land for the 
production of agricultural crops. Hence, people in this area faced 
economic difficulties. To overcome this issue, the people practiced 
sericulture to produce silk as an alternative income source in the 
loft spaces of gassho-style houses [10]. Moreover, according to the 
historical records, the people produced saltpeter in the underfloor 
space of gassho-style houses using ammonia-containing wastes 
such as plant residues (e.g., the stems of Japanese millet, 
buckwheat, and tobacco), silkworm feces from sericulture, and 
human urine as primary substrates [10]. These substrates were 
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mixed with soil under the floor of the gassho-style house to create a 
niter-bed, and nitrate (NO :i ) was produced via a microbial 
nitrification reaction. The resulting nitrate contained in the 
niter-bed soil was further extracted with water and purified by 
adding wood ash to remove unnecessary minerals (e.g., Ca, Mg, 
K, and Na). Purified saltpeter was used as the raw material for 
gunpowder production, and this gunpowder also became an 
important source of income for people who lived in this area 
[9, 1 0] . Thus, saltpeter production via the niter-bed process was a 
cornerstone of the lifestyle in this area. 

Although saltpeter production in the villages of Shirakawa-go 
and Gokayama ended more than 100 years ago, the relict niter- 
beds are now conserved in the underfloor space of gassho-style 
houses, where they are isolated from rainfall and other destabi- 
lizing factors [10]. Surprisingly, bulk soils collected from relict 
niter-beds in Shirakawa-go village retained nitrification activity 
[13]. This finding implies that nitrifying microbes have persisted in 
the relict niter-bed soils; however, there is no information available 
regarding the microbial community structure responsible for 
nitrification in the relict niter-bed ecosystem. In this study, we 
investigated the microbial community structure of nine soil 
samples obtained from the underfloor areas of three gassho-style 
houses by using 454 pyrotag analysis targeting the 16S rRNA gene 
and the ammonium monooxygenase subunit A gene (amoA) of 
ammonia-oxidizing bacteria (AOB) and ammonia-oxidizing ar- 
chaea (AO A). In addition, a previously reported amoA pyrotag 
dataset [14] was used to elucidate the ecological niche of AOB and 
AOA in soil environments. 

Materials and Methods 

Soil sample 

A total of nine soil samples were collected from three relict niter- 
beds under the floors of gassho-style houses located in Shirakawa 
Village, referred to as OVA, OVB, and OVC. We collected five, 
two, and two soil samples from the OVA, OVB, and OVC houses; 
these samples are referred to as OVA1-5, OVB2 and OVB3, and 
OVC1 and OVC2, respectively. Each 500-1,000 g of soil sample 
was collected from relict niter-beds (for OVA and OVC, a depth 
of ca. 10 cm from the surface; for OVB, a depth of ca. 20 cm from 
the surface) by permission from the education board of Shirakawa- 
Village, Gifu, Japan. Detailed physicochemical parameters, 
including pH, organic carbon content, and ion contents, are 
summarized in Table SI [13,15-17]. Pebbles were removed from 
the samples using a 30-mesh filter, and the samples were subjected 
to DNA extraction. 
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DNA extraction, PCR, and pyrosequencing 

DNA was extracted from the soil samples using the FastDNA 
SPIN Kit for Soil (MP Biomedicals, Carlsbad, CA, USA) 
according to the manufacturer's instructions. Approximately 
50 ng of template DNA was used for PCR amplification in a 
reaction volume of 50 |J,L. The 16S rRNA, bacterial amoA, and 
archaeal amoA genes were amplified with the forward/reverse 
primer sets Univ519F/Univ926R [18], amoA-lF/amoA-2R [19], 
and Arch-amoAF/ Arch-amoAR [20], respectively. The forward 
and reverse primers were barcoded with fusion A (FA; 5'- 
CGTATCGCCTCCCTCGCGCCATCAG-3') and fusion B (FB; 
5'-CTATGCGCCTTGCCAGCCCGCTCAG-3') adaptor se- 
quences, respectively. The primers also included Multiplex 
Identifier (MID) sequences. Detailed information regarding the 
primer sequences is provided in Table S2. These genes were 
amplified with PrimeSTAR HS DNA Polymerase (TaKaRa, Otsu, 
Japan) using the following thermocycling conditions: the appro- 
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priate number of cycles (25 cycles for 16S rRNA, 40 cycles for 
bacterial amoA, and 35 cycles for archaeal amoA genes) of 
denaturation at 98°C for 10 sec, annealing at 55°C for 15 sec, and 
extension at 72°C for 1 min. PCR amplicons were purified using a 
Wizard SV Gel and PCR Clean-Up System (Promega, Fitchburg, 
WI, USA) and pooled for subsequent 454 pyrotag analysis. 
Pyrosequencing was performed using the GS-FLX Titanium 
platform (Roche/454 Life Sciences, Branford, CT, USA) at the 
FASMAC Co., Ltd. (Atugi, Japan). Sequencing was performed 
using 1/2-region gaskets. For the 16S rRNA gene, we used 
pyrotag reads obtained from both fusion sides (A and B) because 
the 454 platform used in this study generated full-length amplicons 
(approximately 370 bp). For bacterial and archaeal amoA genes, 
we used pyrotag reads obtained from only the fusion A side 
because the length of each sequenced read did not cover the entire 
region of the amplicon. 

Operational taxonomic unit (OTU) analysis for 16S rRNA 
gene pyrotags 

Raw 16S rRNA gene pyrotag sequences were qualitatively 
screened and trimmed with QIIME 1.7.0 [21], maintaining the 
sequence length at ^300 nt. The sequences had an average 
quality score of &25. A denoising function was applied to the 
filtered sequence data [22]. The denoised sequence data were 
grouped into OTUs with the UCLUST algorithm [23]. The 
sequences exhibiting >97% identities were grouped into one 
OTU. The representative sequences of each OTU were aligned 
with PyNAST [24], and chimeric sequences were detected and 



removed using ChimeraSlayer [25]. The taxonomic assignment of 
each OTU was performed with the RDP classifier trained by a 
dataset from Greengenes at a minimum confidence level of 0.8 
[26]. The sequence data were compared with the 16S rRNA 
database of RDP release 10 [27] using the BLAST+ program [28]. 

OTU analysis for amoA pyrotags 

The procedures for quality screening, trimming, denoising, and 
OTU selection were same as those described above for 16S rRNA 
gene pyrotag analysis, except that sequence identity (>85%) was 
determined according to the threshold of Pester et al. [14]. The 
filtered sequences were manually checked after aligning with 
MEGA program [29]. Chimeric sequences were detected and 
removed using the UCHIME algorithm [30] on the FunGene 
platform [27]. Sequence data were compared with the GenBank 
database using the BLAST+ program [28] . 

Statistical and phylogenetic analysis 

After taxonomy assignments, a-diversity indices (e.g., Chaol, 
Shannon, and Simpson) of the 16S rRNA and amoA genes pyrotag 
libraries were calculated by QIIME [21] and Estimates (version 
9.1.0) [31] software, respectively. The coverage values were 
calculated using equation [1- (n/N)], where n is the number of 
OTUs in a single read (singleton) and N is the total number of 
reads analyzed [32]. The weighted UniFrac distances were 
calculated and used for principal coordinate analysis (PCoA) 
[33]. Phylogenetic trees for the bacterial and archaeal amoA 
libraries for Fast UniFrac analysis were constructed with the 



Group S3: 

GrOUP S2: N0 3 : 2.0-4.6 mg/g 




Figure 1. Principal coordinate analysis based on the abundances of 16S rRNA gene OTUs (weighted UniFrac) and separating soils 
according to their ion contents. For this analysis, observed 16S rRNA gene OTUs were normalized to 7,770 reads per soil. 
doi:10.1371/journal.pone.0104752.g001 



PLOS ONE | www.plosone.org 



3 



August 2014 | Volume 9 | Issue 8 | e104752 



Microbial Community of Relict Niter-Bed in Shirakawa-go 



OTU ID 



OVA 



OVB 



OVC 



• >95 

0>75 
□ >50 



% support 
% support 
% support 




Nitrosospira sp. Nsp5, AY123834 
Nitrosospira sp. Nsp12, AY123823 
NB_Bamo_05 
NB_Bamo_09 

Nitrosospira sp. NpAV, AF032438 
Nitrosospira sp. CT2F, AY1 89143 

1 Nitrosospira tenuis Nv1, AY 123824 

Nitrosospira sp. LT2MFa, AY189145 
NB_Bamo_10 
NB_Bamo_01 

Nitrosospira briensis Nsp10, AY123821 
NB_Bamo_13 
AR_Bamo_01 
Nitrosospira sp. Nsp17, AY123825 
i Nitrosospira sp. Nvo, AY123826 
Jp- Nitrosospira tenuis Nv-12, U76552 
L- NB_Bamo_03 

Nitrosospira multiformis ATCC25196, CP000103 
■AA_Bamo_01 
CR_Bamo_01 

|amoA 20 Haldell kultur, AY795805 
■amoA 27 Sali kultur, AY795811 
Bamo_04 
Nitrosomonas sp. NL7, AY958704 
Nitrosomonas oligotropha, AF272406 
NB_Bamo_08 
Nitrosomonas sp. Nm84, AY123818 
NB Bamo 12 



8 
2' 



NB 



NB 



NB 



NB 



NB 



NB 



Bamo. 
Bamo 
Bamo. 
Bamo. 
Bamo. 
Bamo 



_05O O 
>_09 O 
10 

01 • • 

J3 O 
_03Q O 



o o o o 
o 

o 



o o 




Nitrosomonas sp. JL21, AF327919 

Activated sludge, FV4B4XA02HD703* 
Nitrosomonas ureae, AF272403 
Activated sludge, FV4B4XA02IBDXB* 
r NB_Bamo_11 

* Nitrosomonas marina, AF272405 
— Nitrosomonas aestuarii, AF272400 
-Nitrosomonas marina isolate SF AOB H07, HM345622 



r Nitrosomonas sp. LT-4, JN367456 
NB_Bamo_07 

L- Nitrosomonas europaea ATCC 19178, JN099309 
_r Nitrosomonas sp. GH22, AF327917 
• Nitrosomonas eutropha C91 , CP000450 
NB_Bamo_06 

■ Nitrosomonas communis, AF272399 



AR Bamo 02 



■ Nitrosomonas sp. Nm33, AF272408 
■ Nitrosomonas sp. Nm58, AY123820 
-NB_Bamo_02 

• Nitrosomonas nitrosa, AF272404 



NB 



_Bamo 
_Bamo. 
5-" NB Bamo. 
_Bamo. 
Bamo. 
Bamo 



NB 



NB 



NB 



NB 



08 
12 
11 

'_07O 
06 

02 



o 



o 
o 
o 
o 
o o 
o o o 



i Nitrosococcus oceani Nc10, AJ298699 
' Nitrosococcus mobilis Nc2, AJ298701 



Nitrosococcus 



>50%; 



5-50%; 



0.5-5%; O 0.01-0.5% 



Figure 2. Distance matrix tree of bacterial amoA gene sequences retrieved from niter-bed soils based on the maximum likelihood 
method. Boldface indicates the sequences obtained in this study. The Nitrosococcus amoA gene sequences (AJ298699 and AJ298701) were used as 
outgroups. The bar indicates 10% base substitution. Branching points supported probabilities >95%, >75%, and >50% by bootstrap analyses (based 
on 1,000 replicates, estimated using the maximum likelihood method) are indicated by solid circle, open circles, and open square, respectively. 
doi:1 0.1 371 /journal.pone.01 04752.g002 



MEGA program [29] based on the neighbor-joining algorithm 
[34] . PCoA and jackknife clustering analyses were performed on 
weighted and normalized data with QIIME and the Fast UniFrac 
program [33], respectively. To identify OTUs that exhibited 
significant differences in abundance between different soil samples, 
Welch's i-test (confidence interval method: Welch's inverted, \p< 
0.05] for two groups [definition of each group was described in the 
Results and Discussion]) having possibly unequal variances was 
performed with STAMP software [35,36]. 

The distance matrix tree of the bacterial and archaeal amoA 
genes, which is based on the maximum likelihood method, was 
constructed with the MEGA program [29] . The topology of the 
trees was estimated by 1,000 bootstrap replicates [37]. Previously 
reported amoA pyrotag sequences retrieved from soil environ- 



ments in Austria, Costa Rica, Greenland, and Namibia were 
downloaded from the Sequence Read Archive (SRA) at NCBI 
(accession no. SRA047303) and used for the statistical and 
phylogenetic comparisons [14]. Subclusters within the "Ca. 
Nitrososphaera" and "Ca. Nitrosopumilus" were defined accord- 
ing to Pester et al. [14]. 

Nucleotide sequence accession numbers 

The pyrosequence data obtained in this study have been 
deposited under DDBJ/EMBL/GenBank accession no. 
DRA001202. 
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Figure 3. Principal coordinate analysis based on the abundances of (A) bacterial and (B) archaeal amoA gene OTUs (weighted 
UniFrac) and separating soils according to soil pH and organic carbon contents. For this analysis, amoA pyrotag libraries reported by 
Pester et al. [14] were used for the comparison, and bacterial and archaeal amoA OTUs observed here were normalized to 2,320 and 2,130 reads per 
soil, respectively. 

doi:1 0.1 371 /journal.pone.01 04752.g003 



Results and Discussion 

16S rRNA gene pyrotag analysis 

In this study, nine soil samples from three gassho-style houses, 
referred to as OVA (samples OVA1-5), OVB (samples OVB2 and 
OVB3), and OVC (samples OVC1 and OVC2), were used for 
microbial community analysis using 454 pyrosequencing. We 
retrieved a total of 173,016 16S rRNA gene pyrotags and found 
3,146 OTUs with the criterion of >97% sequence identity 
(Table 1). Although the rarefaction curves showed that the 
accumulation of the 16S rRNA gene OTUs was insufficient to 
achieve the goal of saturation (Fig. SI A), the coverage values were 
calculated to be 96.99-98.59%, suggesting that the OTUs 
retrieved here were sufficient for estimating the microbial diversity 
in the relict niter-bed ecosystems. According to the Chaol 
nonparametric estimators, the soil samples contained approxi- 



mately 1.38-1.63-fold more OTUs than what was detected in our 
analyses. 

The OTUs belonging to the phyla Proteobacteria (especially the 
classes Alphaproteobacteria and Gammaproteobacteria), Actinobac- 
teria, Bacteroidetes, Chloroflexi, Firmicutes, Gemmatimonadetes, 
and Planctomycetes predominated and were detected in all samples 
(Table S3 and S4). Such microbial taxa are frequently found in soil 
environments [38]. In each sample, only a few number of the 16S 
rRNA gene pyrotags (<0.8% of the total populations) were 
associated with previously known ammonia-oxidizing microor- 
ganisms at the genus level: for AOB, three OTUs (nos. 1267, 2435, 
and 2701) of the Nitrosococcus , one OTU (no. 2080) of the 
Nilrosomonas, and one OTU (no. 673) of the Nitrosospira; for 
AOA, four OTUs (nos. 1207, 1897, 1946, and 3119) of the "Ca. 
Nitrososphaera" (Table S5). Given that soil samples from the relict 
niter-beds have retained microbial nitrification activity [13], these 
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Figure 4. Distance matrix tree of archaeal amoA gene sequences retrieved from niter-bed soils based on the maximum likelihood 
method. Boldface indicates the sequences obtained in this study. The amoA gene sequence of the "Ca. Nitrosocaldus yellowstonii" strain 
HL72 (EU239961) was used as outgroup. The bar indicates 10% base substitution. Branching points supported probabilities >95%, >75%, and >50% 
by bootstrap analyses (based on 1,000 replicates, estimated using the maximum likelihood method) are indicated by solid circle, open circles, and 
open square, respectively. 
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nitrifying microbes may contribute to the nitrification at low 
population densities; alternatively, there may have been uniden- 
tified nitrifying microbes in the relict niter-bed soils. Although the 
relative abundances of nitrifying microbes were in good agreement 
with the microbial compositions of several types of soil ecosystems 
[38], we could not fully elucidate the biodiversity and ecological 
niche of ammonia-oxidizing microbes in the relict niter-bed 
ecosystems based on the results of 16S rRNA gene pyrotag analysis 
(e.g., there are no significant relationship between the ratio of the 
abundance ratio of AOB to AOA and clustering of the 
communities in PCoA and jackknife analysis described below) 
(Table S5). These observations emphasize the need to implement 
further analysis targeting amoA gene that enable specific detection 
of ammonia-oxidizing microbes. 

The similarity of microbial compositions among the nine soil 
samples was compared on the basis of 16S rRNA gene pyrotag 
sequences. According to the PCoA and jackknife analyses, the 
microbial community structures of the relict niter-bed soil were 
classified into three groups (Fig. 1 and S2A). Group SI, consisting 
of OVA1 and OVA2, was distinct from the other two groups. 
Based on previously reported physicochemical characterization 
[15-17], relatively high contents of N0 3 ~, Cl~, K + , and Na + were 
found in OVA1 and OVA2 compared with the samples clustered 
into groups S2 and S3 (Table SI). Because the soils in group SI 
were not subjected to the extraction step of saltpeter production, 
the ion contents were higher than those in the other relict niter- 
bed soil samples [17]. In addition, the N0 3 and CI contents 
were relatively low in the soils clustered in group S3 compared 
with those in group S2. Based on Welch's t-te&t, we found a total of 
21 OTUs that were significantly different in relative abundances 
among the groups and were associated with dominant taxa (e.g., 
Alphaproteobacteria, Gammaproteobacteria, Actinobacteria, Bacter- 
oidetes, Chloroflexi, Firmicutes, and Gemmatimonadetes) (Fig. S3). 
Among them, Gammaproteobacteria-type OTU (NB_16S_507) 
and Chloroflexi-type OTU (NB_16S_1415) were involved in the 
differentiation of groups SI, S2, and S3. The samples in group S2 
were further separated according to their moisture contents 
(groups S2-1 and S2-2), and the relative abundances of six OTUs 
(NB_16S_1093, 1144, 1409, 1415, 2046, and 2058) were 
significantly affected by the moisture contents (Fig. S3). Several 
types of ion species and moisture contents have been recognized as 
potential factors in the determination of soil microbial constituents 
[39-42]. Thus, the NO :i ", CI", K + , Na + , and moisture contents 
may affect the microbial community structures of relict niter-bed 
ecosystems. 

Bacterial amok gene pyrotag analysis 

As described above, the proportions of 16S rRNA gene pyrotags 
associated with AOB and AOA were very low in the relict niter- 
bed soils. To further elucidate the microbial constituents 
responsible for ammonia oxidation, we performed pyrotag analysis 
targeting bacterial and archaeal amoA genes. As a result, we 
obtained a total of 59,404 bacterial amoA pyrotags and found 13 
OTUs in relict niter-bed soils (Table 1). However, the bacterial 
amoA gene fragment could not be amplified from OVA1 and 
OVA2 using the primer set described in this study, suggesting that 
there were either low levels of AOB in the samples or there were 
some AOB with significantly different types of amoA genes. 
Instead, the archaeal amoA gene fragments were successfully 
amplified from OVA1 and OVA2 samples, suggesting that AOA 
may be a major player in ammonia oxidation in these soils (details 
described below). Although the rarefaction curves for the bacterial 
amoA libraries shorted to reach a plateau (except for OVC1 and 
OVC2), the coverage values were calculated to be 99.96-100%, 



indicating that the OTUs retrieved here were sufficient for 
estimating the biodiversity of AOB in seven relict niter-bed soils 
(Table 1 and Fig. SIB). 

The most abundant AOB belonged to the genus Nilrosospira (6 
OTUs, 68.13-99.96% of the total reads in each sample), and 
Nitrosospira briensis-type OTU (NB_Bamo_01) predominated in 
all samples (except for OVA1 and OVA2, in which there were no 
bacterial amoA gene amplification) (Table S6). This result is in 
good accordance with 16S rRNA gene pyrotag analysis that found 
a Nitrosospira briensis-type OTU (no. 673) across the samples. An 
OTU (NB_Bamo_01) related to the Nitrosospira briensis strain 
NsplO (95.5% amoA gene sequence identity) was detected as a 
major AOB population in all niter-bed soils (Fig. 2 and Table S6). 
Nitrosospira briensis are enriched in and have been isolated from 
soil samples taken from sandy and agriculturally uncultivated areas 
of Europe [43] . To the best of our knowledge, only two reports 
have described the clear appearance of Nitrosospira briensis 
populations: water columns of the Chesapeake Bay [44] and soil 
samples of the Cascade Mountains in Oregon [45] in the U.S. 
Interestingly, Bollmann et al. [46] reported that under ammonia- 
starved conditions, the cells of Nitrosospira briensis had the ability 
to immediately recover ammonia oxidation activity upon exposure 
to ammonia-rich conditions. Although the mechanism behind this 
phenomenon is unclear, Nitrosospira briensis-type AOB may 
reserve their activity in ammonia-limited environments, such as 
the relict niter-bed soils. Within the genus Nitrosospira, an OTU 
(NB_Bamo_03) that is closely related to the Nitrosospira tenuis 
strain Nv-12 (96.2% sequence identity) was also detected in all 
bacterial amoA gene pytotag libraries; however, its ecophysiolog- 
ical traits (other than ammonia oxidation) remain unclear. 

In addition to the Nitrosospira populations, six Nitrosomonas- 
type OTUs were found in seven niter-bed soils (Fig. 2 and Table 
S6). In contrast to the Nitrosospira-type OTUs, 16S pyrotags 
related to Nitrosomonas were found in OVA3 and OVA4; 
however, the proportions were quite low (0.005-0.006% of total 
16S rRNA gene pyrotags). This finding may also support the 
utility of amoA -based pyrotag analysis for relict niter-bed soils. An 
OTU (NB_Bamo_02) associated with Nitrosomonas nitrosa 
(87.2% sequence identity) was observed with significant abun- 
dance in OVC1 and OVC2 (5.6% and 31.8%, respectively) (Table 
S6). According to previous reports, Nitrosomonas nitrosa-type 
AOB have frequently been found in estuary soils [47], nitrification 
bioreactors [48,49], and freshwater ecosystems [50]. In addition to 
these taxa, an OTU (NB_Bamo_04) related to the unidentified 
amoA gene cluster consisting of clones retrieved from potash- 
polluted marsh soil [51] was found in three relict niter-bed soils 
(Fig. 2 and Table S6). This cluster has a significantly low identity 
with previously reported bacterial amoA gene sequences (e.g., 
NB_BamoA_04 had 77.4% identity with Nitrosospira multiformis 
strain C-71, accession no. X90822), suggesting that phylogenet- 
ically novel AOB may play a role in ammonia oxidation in some 
niter-bed ecosystems. 

The similarity of AOB compositions in relict niter-bed soils was 
compared with previously reported amoA gene-based microbial 
community analysis [14]. PCoA and jackknife analyses showed 
that soil bacterial amoA libraries were separated into two groups 
according to their pH: seven niter-bed soils were clustered into 
group Bl with two Austrian arable and riparian forest soils (pH 
7.05-8.35), and the remaining two soils, from an Austrian spruce 
forest and a Costa-Rican rain forest, were clustered into group B2 
(pH 4.36-4.63) (Fig. 3A). pH is generally to be a key parameter in 
defining the abundance of AOB populations in soil ecosystems 
[52,53] and bioreactors [54,55]. The niter-bed soils in group Bl 
were further divided into two subgroups (Fig. S2B and S3A); 
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however, there was no clear explanation for the differentiation 
between OVC2 and the other soils based on the physicochemical 
parameters analyzed herein. In short, soil pH is one of the 
important factors affecting the AOB populations in relict niter-bed 
ecosystems. 

Archaeal amoA gene pyrotag analysis 

A total of 33,087 archaeal amoA gene sequences were retrieved 
from nine relict niter-bed soils, and we found 16 OTUs with the 
criterion of >85% sequence identity. The archaeal amoA diversity 
was almost fully covered according to the rarefaction curves, 
coverage values, and Chaol diversity index (Table 1 and Fig. 
SIC). 

Members of the genus "Ca. Nitrososphaera" were detected as 
major AOA populations (12 OTUs, 99.23-100% of the total reads 
in each sample), whereas those of "Ca. Nitrosopumilus" were 
absent or few (<0.28% of the total reads in each sample) in all 
samples tested (Fig. 4 and Table S6). Based on the 16S rRNA gene 
pyrotag analysis, only a few OTUs (accounting for 0.16-1.1% of 
the total prokaryotic population) associated with "Ca. Nitroso- 
sphaera" were identified as AOA populations in relict niter-bed 
soils (Table S5). Thus, the results from amo^l-based pyrotag 
analysis expanded the biodiversity of AOA in relict niter-bed soils. 
The most abundant OTU in OVA1 and OVA2 soils was 
NB_Aamo_02, which was categorized into subcluster 4.1 of 
"Ca. Nitrososphaera" and was distandy related to previously 
characterized AOA (77.8% sequence identity with "Ca. Nitroso- 
sphaera viennensis" strain EN76, accession no. FR773159), 
whereas the most abundant OTU in the other seven samples 
(NB_Aamo_01) was closely related to "Ca. Nitrososphaera 
gargensis" strain Ga9.2 (90.6% sequence identity, accession no. 
CP002408) and categorized into subcluster 1.1. This marked 
difference in a major AOA population within OVA soils was 
clearly associated with the maximum nitrification rate of bulk soil 
(Table SI) [13]. This finding also implies that the nitrification 
activity and/or affinity for ammonia may differ between subclus- 
ters 1.1 and 1.4. The amoA fragments associated with "Ca. 
Nitrososphaera" subcluster 1.4 have been observed in diverse 
environments, including soil [20,56,57], salt marsh sediment [58], 
and marine sediment [20,59]. In addition to the major AOA, 
several subclusters of "Ca. Nitrososphaera" were detected as 
minor AOA populations. For example, an OTU (NB_Aamo_03) 
assigned to subcluster 9 was commonly observed in the relict niter- 
bed soils (Fig. 4). Although information regarding the ecological 
traits of "Ca. Nitrososphaera" subcluster 9 organisms is limited, 
these organisms were widely distributed in natural environments, 
such as marine sediment [60] and soil [20]. Further studies using 
metagenomic and single-cell genomic approaches are necessary to 
elucidate the ecological functions of unidentified AOA in relict 
niter-bed ecosystems. 

PCoA and jackknife analyses clearly showed that archaeal amoA 
pyrotag libraries were separated into three clusters according to 
soil pH and organic carbon concentration (Fig. 3B and S2C). The 
niter-bed soils formed group Al with an Austrian riparian forest 
soil. Group Al soils had a weak alkaline pH (7.05-8.35) and 
relatively high organic carbon content (3.7-23.2%). In addition, 
relict niter-bed soils were further separated into two subgroups 
according to their N0 3 ~, Cl~, K + , and Na + contents in addition 
to the 16S rRNA gene-based comparison (Fig. S4B). The samples 
clustered into group A2 were mainly composed of Namibian soils, 
which had relatively low organic carbon contents (0.29-1.61%) 
and a weakly acidic pH (5.7-6.8). The samples categorized into 
group A3 (soil samples from a Costa-Rican rain forest, a Costa- 
Rican arable forest, and Greenland tundra) were defined as having 



a relatively high organic carbon content (1.66-7.51%) and an 
acidic soil pH (4.63-5.09). Interestingly, there was a significant 
correlation between the Chao 1 diversity index and organic carbon 
content, with higher richness at lower levels of organic carbon (Fig. 
S5A). Moreover, a relatively weak correlation between the Chaol 
index and soil pH was also observed, with the highest diversity at a 
pH of approximately 6 (Fig. S5B). These results were in good 
agreement with the previously reported amoA gene-based pyrotag 
analysis [14], and they provided further evidence of the 
importance of soil pH and organic carbon for defining the 
ecological niche of the AOA community in soil ecosystems, 
including relict niter-beds. 

Conclusions 

16S rRNA and amoA gene-based pyrotag analyses revealed that 
microbes responsible for ammonia oxidation (both AOB and 
AOA) still persist in relict niter-bed soils, although saltpeter 
production ceased more than 100 years ago. Our study provided 
the first insights into the microbial diversity in relict niter-bed 
ecosystems. Our main conclusions are as follows, (i) Based on 16S 
rRNA gene pyrotag analysis, members of the phyla Proteobacteria, 
Actinobacteria, Bacteroidetes, Chloroflexi, Firmicutes, Gemmatimo- 
nadetes, and Planctomycetes were major microbial constituents in 
the niter-bed soil samples. The N0 3 , CI , K + , Na + , and moisture 
contents were potential factors for the microbial compositions of 
the relict niter-bed soils, (ii) Members of Nitrosospira and "Ca. 
Nitrososphaera" were the major AOB and AOA populations, 
respectively, in the niter-bed soils. Soil pH, organic carbon 
content, and N0 3 ~, Cl~, K + , and Na + contents were important 
factors for the ecological niche of AOB and AOA in the relict 
niter-bed soils. 

Supporting Information 

Figure SI Rarefaction curves of (A) 16S rRNA, (B) 
bacterial amoA, and (C) archaeal amoA gene sequences 
of nine niter-bed soils. 

(TIF) 

Figure S2 Jackknife clustering of (A) 16S rRNA, (B) 
bacterial amoA, and (C) archaeal amoA gene pyrotag 
libraries from the relict niter-bed soils (weighted 
UniFrac, normalized to 7,770, 2,320, and 2,130 reads 
per soil, respectively). 
(TIF) 

Figure S3 Extended error bar plot showing the OTUs 
that have significantly (p<0.05) different abundances in 
the group. Groups were defined on the basis of principle 
coordinate analysis (Fig. 1) and Jackknife clustering of 16S rRNA 
gene pyrotag analysis (Fig. S2A). The OTUs overrepresented in 
upper group have a positive difference between relative abun- 
dances, and those overrepresented in bottom group have a 
negative difference between relative abundances. Only OTUs with 
mean proportions more than 0.01% in the total populations were 
shown in the plot. 
(TIF) 

Figure S4 Principal coordinate analysis based on the 
abundances of (A) bacterial and (B) archaeal amoA gene 
OTUs (weighted UniFrac) only for niter-bed soils. For this 
analysis, bacterial and archaeal amoA OTUs observed here were 
normalized to 2,320 and 2,130 reads per soil, respectively. 
(TIF) 
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Figure S5 Correlation analysis of organic carbon con- 
tents and soil pH to Chaol diversity indexes. Previously 
reported archaeal amoA pyrotag data from various soil environ- 
ments [14] were combined with those from niter-bed soils. 
(TIF) 

Table SI Physicochemical parameters of niter-bed soil 
samples. 

(PDF) 

Table S2 Primer sequences for pyrotag analysis. 

(PDF) 

Table S3 The 1 6S rRNA gene pyrotag libraries of relict 
niter-bed soil. 

(PDF) 

Table S4 The OTU table of 16S rRNA gene pyrotag 
libraries. 

(PDF) 
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